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1. Introduction 

A major issue of Lattice Gauge Theory has been the reduction of effects which are due to the 
finite lattice spacing a, in order to better approach the elusive continuum limit. Over the years, 
many efforts have been made for ^(a^) improvement in lattice observables, which in many cases 
is automatic by virtue of symmetries of the fermion action. According to Symanzik's program |jl|], 
one can improve the action by a judicious addition of irrelevant operators. Also, in the twisted mass 
formulation of QCD at maximal twist, certain observables are ^(a^) improved, by symmetry 
considerations. The first 1-loop perturbative computation of 0'{a^) effects was recently performed 
by our group [^. This regards the evaluation of the fermion propagator and Green's functions for 
ultralocal bilinears of the form = ^TXa^, using the Wilson/clover fermions and Symanzik 
improved gluons. Extending a calculation up to brings in new difficulties, compared to 

lower order in a; for instance, there appear new types of singularities. The procedure to address 
this issue is extensively described in Ref. [^]. The G{cP-) terms of such perturbative computations 
are of great utility since they can be subtracted from non-perturbative estimates to minimize their 
lattice artifacts. 

The generalized parton distributions (GPDs) of the nucleon determine non-forward matrix el- 
ements of gauge invariant light cone operators; the moments of such operators can be evaluated 
on the lattice. GPDs also give information on interesting quantities such as the quark orbital an- 
gular momentum contribution to proton spin [Q]. Moreover, the generalized parton distributions 
can be measured in high energy scattering experiments, for instance the deeply virtual Compton 
scattering of virtual photons off a nucleon. Thus, by computing moments of GPDs in lattice QCD 
one can explore many aspects of the nucleon structure. The operators that are related to the GPDs 
must be renormalized, before one compares results from simulations to physical, experimentally 
measurable quantities. 

In this work, we investigate the perturbative renormalization of the fermion propagator, local 
and twist-2 fermion operators. We compute all matrix elements for the amputated Green's functions 
of the inverse fermion propagator, ultralocal bilinears and twist-2 one-derivative operators, defined 
as ^^Tj^ D v}T"*P (symmetrized and traceless). Although our expressions for all these matrix 
elements are as general as possible (the dependence on almost all parameters is shown explicitly), 
these are extremely lengthy and complicated to be presented here. Thus, we show our results only 
for the renormalization constants and for particular choices of the various parameters. In Section ^ 
we briefly describe the procedure for these computations and in Sections ^ - ^ we present our results 
for the renormalization functions. We also compare with non-perturbative estimates provided by 
the ETM Collaboration. Due to space limitations the renormalization of the ultralocal operators is 
not discussed here, but will appear in a longer write-up. 



2. Description of the calculation 



In the framework of this computation we employ a fermion action which includes the clover 
term, and also an additional mass parameter, \i 
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A summation over the Dirac indices V, p is implied. The advantage of the above action is that it 
combines different actions that are used widely: Wilson (csw — )• 0, /i — )• 0), clover (/i 0) and 
twisted mass (csw 0). It is important to mention that for the case of twisted mass fermions, 
Eq. (2J^) corresponds to the twisted mass action in the so called twisted basis % {% = e^'T^s^^xj/) 



The appropriate form of this action in physical basis has a rotation of the Wilson term, r — > 
—iry^T^. This is crucial because the definition of the Z-factors in physical and twisted bases 
are different, for instance Z'^-^ = z'y'^"^'^ and Zy^^^ = z'^''^"^"^ , for the axial and vector currents with 
or without a derivative. The results presented here are all in the physical basis. 

For the gluon part we use a family of Symanzik improved gluons and more precisely 10 dif- 
ferent sets of Symanzik parameters that are widely used in simulations; these are shown in Table 1 
of Ref. |j3|]. In the same reference, the reader can find all details for the calculation procedure and 
the evaluation of the superficially divergent integrals. The substantial difference with our previous 
computations is the presence of masses, m and jU, in the fermion propagators. The dependence 
of superficially divergent integrals on p,m, jj. is nontrivial and it is a rather complicated task to 
extract this dependence explicitly. In fact, there are 11 "strongly" IR divergent integrals (they are 
convergent only at D > 6 dimensions) and a few hundred of divergent integrals at D = 4 — 2e 
dimensions. 

The 1-loop diagrams that enter our three calculations are shown in Fig. 1. Diagrams la and 
lb regard the correction to the inverse fermion propagator, diagram 2 corrects the local bilinears 
and finally diagrams 3a-3d contribute to the calculation of the twist-2 bilinears. The appropriate 
operator employed in each computation is represented by a cross in diagrams 2, 3a-3d. 




la 



lb 




3a 3b 3c 3d 

Figure 1: One-loop diagrams contributing to the correction of the amputated Green's functions 
of the propagator (la, lb), local bilinears (2) and twist-2 operators (3a-3d). A wavy (solid) hne 
represents gluons (fermions). A cross denotes an insertion of the operator under study. 

3. Renormalization of the fermion field to i&(a^) 



For all the renormalization factors used here we employ the RI'-MOM scheme. The fermion 
field renormalization constant, Zq, can be obtained using alternative renormalization conditions that 
differ in their lattice artifacts. The most widely used conditions are the following two 



1 

4^' 



1 

4^' 
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(3.1) 
(3.2) 
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with the trace taken over the Dirac indices. The S^\^^{p) is the amputated Green's function for 
the inverse propagator up to 1-Ioop; it was computed from diagrams la and lb of Fig. 1. Indeed, 
the constants and differ only in their lattice artifacts, since the ^(a^) expansion of 5'tree(/') 

is simply S^on'L{p) + ia^ LpP^ Jp/i^P^)- 

In the most general expression for Zq, we show explicitly the dependence on the action pa- 
rameter csw, the coupling constant g, the number of colors Nc, the gauge fixing parameter A, the 
masses m,ix, the lattice spacing a and the external momentum p. On the contrary, we cannot 
express Zq in a closed form as a function of the Symanzik parameters, cq, c\ , C2, cy, we have com- 
puted Zq for each of the 10 sets of Symanzik coefficients shown in Table 1 of Ref. [^]. Next we 
provide our result for Zq using the condition in Eq. (3J), for the Landau gauge, csw = 0, m = 
and tree-level Symanzik improved gluons 
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(3.3) 



In the above expression ;?4 = LvPv" = i^c ~ l)/(2A^c), and is the renormalization scale. 



No ff{a^) terms appear in Eq. (3.3) since we set m equal to zero. Although the mass dependence 



shown in Eq. (3.3) is very complicated, the numerical result for the renormalization constant Zq 
depends mildly on /i (for the values used in numerical simulations by ETMC, 0.003 < jX < 0.01). 

A very important issue is that the ^(a^) terms depend not only on p^, but also on the direction 
of the momentum p, as manifested by the presence of pA. As a consequence, different renormal- 
ization prescriptions, involving the same renormalization scale fx but different directions of p, lead 
to a different renormalization constant. The strategy one may follow is to average over the results 
for Zq which are obtained using these different directions. 

In Table 1 we present our 1-loop results for Zq, up to G{a?-\ for the two renormalization 



conditions: Eq. (3J) (Z^) and Eq. ( ^^ ) (Z^). The results are given in the chiral limit for the 
Landau gauge, csw = 0, Nc = 3 and for 3 values of the coupling constant j3 = 2Nc/g^- Instead of 
the bare coupling we have used the tadpole-improved coupling [^, defined as = j {plaq), to 
achieve further improvement. The renormalization scale }X was set to \/a, that is a^p^ = 1. This 
is a well-defined choice, since all possible lattice momenta (permutations of (±1,0,0,0)) have the 



same pA. Thus, Zq is the same for all these permutations. Z^^^" and Zq'^'^" ' are the ff{a^) effects 
of Zq and Zq respectively. The last column gives the non-perturbative estimates of ETMC p]. 



using the renormalization condition of Eq. (3.2) and employing the same pai^ameters as we did. In 



^non-pert 



our 0'{a^) corrections were already subtracted to reduce lattice artifacts. In fact, ETMC 



applied the subtraction procedure to all their data; therefore, the behavior of Z; 



-non— pert 



against the 
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renormalization scale becomes flatter in the energy region where perturbation theory is valid. This 
indicates that the lattice artifacts are suppressed. 



i8 








^q 


rynon—pert 
^q 


3.80 
3.90 
4.05 


0.655255014(4) 
0.663871997(4) 
0.675276827(4) 


-0.039554503(1) 
-0.045705298(1) 
-0.053846057(1) 


0.771056601(4) 
0.782134880(4) 
0.796797308(3) 


0.076247083(1) 
0.072557584(1) 
0.067674424(1) 


0.755(4) 

0.757(3) 
0.777(5) 



Table 1: Perturbative and non-perturbative estimates of for different coupling constants, in the chiral 
limit and using two renormalization conditions, for tree-level Symanzik gluons, csw — Q,Nc = 3. 



The numbers is parenthesis denote the systematic errors coming from our procedures for numerical 
integration over loop momenta. Since we have used two alternative renormalization conditions, it is 
worth comparing their lattice artifacts. For p.^ =4 and the same parameters of Table 1, we isolated 
the contribution of the 6(0^) terms and tabulated them in Table 2. In the same Table we present 
the averaged results, {Zq^^" ■*), {Zg^^" ■*). For the averaging we used all 24 lattice momenta with 
a^p^ = 4: 16 permutations of (ibl,zbl,ibl,±l) and 8 permutations of (±2,0,0,0). One observes 
that, although the condition in Eq. ( ^ ) gives ^{a^) correction at tree level for Z^, the combined 
G(a^) effects in Z^ are much smaller than the ones of Z^, which is exactly 1 at tree level. It is also 
evident that averaging over momenta with different directions helps to reduce the overall ff{a^) 
contributions. 



i3 


^q 


(zf'("')) 


ryB0(a^) 

^q 




3.80 
3.90 
4.05 


0.35679267(9) 
0.26014128(8) 
0.13958988(7) 


0.29105315(4) 
0.23208860(4) 
0.15854326(3) 


0.45876742(9) 
0.41544327(8) 
0.36140590(7) 


0.34204053(4) 
0.30973960(4) 
0.26945127(3) 



Table 2: The ff[a'^) contributions up to 1-loop in the perturbative results for the Zq at /x^ = 4. 



4. ^(a^) corrections to the renormalization of twist-2 fermion bilinears 



In this section we present the computation of the amputated Green's functions for the following 
two twist-2 operators 



5){Vl,V2} _ 1 

^ " 2 



5{Vl,V2} 



1 



1 



-5,,,,Y.^ypDpZ^^> (4.1) 

p 

\5y,y,Y.^ypDpx''^ (4.2) 
^ p 



which are symmetrized and traceless, to avoid mixing with lower dimension operators. We have 
computed, to ff{a^), the matrix elements of these operators for general external indices Vi, V2, 
mass, g, Nc, a, p, csw and A. The final results are available for the 10 sets of Symanzik coefficients 
we have used in the calculation of Z^. 
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The one-derivative operators fall into two different irreducible representations of the hypercu- 
bic group, depending on the choice of the external indices Vi , V2. Hence, we distinguish between 



7Vi=V2 



72 



. 7Vi=V2 



72 
^A 



(4.3) 



The renormalization conditions from which Zy and Za are obtained, is defined in the RI'-MOM 
scheme, as 



7 7 



1„V,V2 



Tr 



^ff \P) ^^ffcont 



Tr 



.V1V2 .V,V2 

ffcont ffcont 



AV1V2 
'^ffcont 



where = 7vi , 75 7vi- The quantity A^^^(p) is our result up to 1-loop for the amputated Green's 
function for each operator. In the renormalization condition there is a choice of using A^^^^^ in- 
stead of A^^^^jj. The final Z-factors of the two choices have different lattice artifacts. Applying the 
renormalization condition and using the Zq that we computed perturbatively with its ^{a^) correc- 
tions, we have the results for Zy and Za- In Table 3 we present these results in the chiral limit, for 
tree-level Symanzik gluons, j8 = 3.9, csw = and Landau gauge. We chose six different values of 



the momentum, thus six renormalization scales. For Z^ we employed Eq. (3J), and we used A^^^^^^ 



in Eq. (p!4|). The extrapolation errors are smaller then the last digit shown in Table 3. 



a p 


7I 


72 


7I 

^A 


72 

^A 


(4,2,2,2) 


1.04129 


1.05738 


1.08270 


1.06309 


(5,2,2,2) 


1.04411 


1.04331 


1.08772 


1.04669 


(6,2,2,2) 


1.04881 


1.02365 


1.09312 


1.02496 


(3,3,3,3) 


1.00343 


1.05427 


1.03778 


1.03607 


(4,3,3,3) 


0.99660 


1.05237 


1.04456 


1.03289 


(5,3,3,3) 


0.99731 


1.04736 


1.05543 


1.02660 



Table 3: Perturbative results for Za and Zy in the chiral limit, for various renormalization scales. 
The action parameters are: tree-level Symanzik gluons. Landau gauge, j3 — 3.9, csw = 0, Nc = 3. 



For the case of the vector operator with Vi / V2 we give the result for Tr A^'^^ (p) ■ A^! 
contributes to Zy, for the Landau gauge, tree-level Symanzik gluons, csw = 0, /X = 



1V2 

f cont 



which 



Tr 



V2 

cont 



= -p 



iPti+pti , Cf8^ 



Vi FV2 



15.045752(l)p 



2 , 4pt,Pv2 



3p2 



8p'v, Hp') 



-0A32l2{3)p'-pi -4.6352(l)p4 -4.0096(1);4, 



3 ■ y .'V, 

353 p4 752^ 29 P^pI.pI, , 179 (4, pi, +Pv, 
720 



'ViPv2 
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' Vi V2 



, , 2 103 2 2 331 4 

^ ^ 360 ' ' 360 ^ 



symmetric terms : Vi O V2 



90 p 
1013 



180 



2 2 

Pv, PV2 



(4.5) 
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It is now interesting to use non-perturbative estimates for Za and Zy, combined with our 0'{a^) 
terms, in order to see if the subtraction procedure works well for these operators. Our collaborators 
in ETMC provided us with data for the physical [Q] corresponding to the six momenta of Table 
3. Note that the momenta can be grouped into two categories, according to their spatial values: 
(2,2,2) or (3,3,3). Non-perturbatively it is observed that these groups have different lattice artifacts 
and a discontinuity appears in the Z-factors with one derivative, as shown in Fig. 2. We also plot 
the subtracted non-perturbative data, which exhibit a smooth behavior of Z^. The errors are too 
small to be visible. 



1.25 



1.2 



<; 

N 1.15 



1.1 



1.05 



1.2 1.4 1.6 l.i 
2 2 

a p 



unsubtracted 
subtracted 



2.2 2.4 



Figure 2: Physical non-perturbative with subtractions of 6{a^) terms. 



Another plot demonstrating the achieved improvement of non-perturbative estimates for the 
renormalization constants using our correction terms, is shown in Ref. [^]. In particular, for the 
local axial renormalization constant it is shown that, by converting the subtracted numbers into MS 
scheme and setting the renormalization scale to 2GeV, one obtains a good plateau when plotting 
Za for each of the momentum 4-vectors. 
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